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There is compelling evidence from inelastic-neutron-scattering and tunneling experiments that the 
heavy-fermion superconductor UPd2Al3 can be understood as a dual system consisting of magnetic 
excitons, arising from crystal-field-split U 4+ levels, coupled to delocalised f-electrons. We have 
computed the superconducting transition temperature and the mass renormalisation arising from a 
(N ■ dual model with maximal spin anisotropy using a strong-coupling approach. We find an instability to 

two possible opposite-spin-pairing states with even- or odd-parity gap functions. Each has a line node 
perpendicular to the c-direction, in agreement with NMR relaxation-rate, specific-heat and thermal- 
D ' conductivity measurements. In addition, both have total spin component 5*2=0, compatible with the 

»' | , observation of a pronounced Knight shift and H C 2 Pauli limiting. For parameter values appropriate 

to UPd2Al3, we determine the dependence of the superconducting transition temperature T c on a 
phenomenological coupling constant g and we investigate the associated mass enhancement and its 
anisotropy. 



PACS numbers: 74.20.Mn, 74.20.Rp 

g 

O I- INTRODUCTION 

O ' 

Heavy-fermion superconductors are complicated materials which possess a fascinating and rich variety of physical 
properties and which have stimulated the creation of an equally diverse range of theories in an effort to understand 
them. These materials are dominated by strong electronic correlations giving rise to magnetic and superconducting 
, instabilities. Accordingly, most theories ignore phononic degrees of freedom. 

One of the most extensively studied models of heavy-fermion superconductivity is a phenomenological model based 
on the exchange of spin fluctuations between heavy quasiparticles. The latter originate in periodic resonant scattering 
of conduction electrons close to the Fermi level as described within the Kondo lattice model of Ce compounds 1 ' 2 . 
' Pairing via spin fluctuations can be considered as a one-component model to the extent that the spin fluctuations 
originate in the system of heavy quasiparticles and any interaction between conduction electrons and more localised 
electrons, or conduction electrons and phonons, is neglected. In spite of its simplicity this model has proved to 
be very useful in our attempts to understand heavy-fermion compounds on the border of magnetic long-range order, 
such as CePd2Si2 3 , Celn3 3 and CeRhlns 4 . Its success probably derives from the fact that the spin-fluctuation-induced 
effective interaction tends to dominate all other channels of interaction when a material is tuned close to the border 
of magnetism by doping or, as in the Ce compounds mentioned above, by applying pressure. The spin-fluctuation 
mechanism in the Kondo lattice may also be appropriate for the Ce-based heavy-fermion superconductors at ambient 
Q . pressure, CeCu2Si2 and CeCoIn 5 4 ' 5 . 

It has become increasingly evident that the Kondo-lattice model is, however, not adequate in the case of heavy- 
quasiparticle formation in some Uranium heavy-fermion compounds 6 - , where 5f electrons are partly localised and 
partly itinerant. The former occupy crystalline-electric-field (CEF) split 5? states. The latter, more itinerant, 5f 
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electrons have a strongly enhanced effective mass due to a coupling to virtual excitations between CEF states. This 
notion of heavy-quasiparticle formation is especially appropriate for UPd2Al3 and it is the starting point of our 
theoretical model for heavy-fermion superconductivity in this compound. Specifically, we consider a novel effective 
pairing mechanism which is based on the virtual exchange of collective CEF excitations known as 'magnetic excitons'. 
They are propagating bosonic modes in contrast to the overdamped modes which give rise to pairing in the itinerant 
spin-fluctuation model. 

The interplay of conduction electrons and CEF excitations is a well studied subject. Initially, research was focused in 
two areas 8 : transport anomalies and superconducting pair breaking or enhancement due to conduction-electron scat- 
tering from dilute CEF-split impurities; and periodic lattices of CEF ions interacting with one another via conduction- 
electron polarisation. In the latter case attention was focused on collective effects within the CEF system due to the 
RKKY interactions mediated by itinerant electrons. Later, White and Fulde 9 showed that the inverse effect, viz. 
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conduction-electron mass enhancement via virtual emission and absorption of magnetic excitons, is also important. 
They demonstrated that this mechanism explains the enhanced effective electron masses in praseodymium metal. 

It is then natural to ask whether the exchange of magnetic excitons between quasiparticles can also mediate 
superconductivity. Until recently, no theoretical work had been carried out to answer this question, mainly for two 
reasons. Firstly, when the effects on superconductivity of paramagnetic impurities with CEF-split energy levels were 
studied in detail 10 , it was found that for rare-earth ions the pair-breaking transition matrix elements usually dominate 
the pair-enhancing matrix elements. Secondly, no good example of a superconducting compound with this dual nature, 
i.e., magnetic excitons arising from localised f electrons coupled to delocalised f electrons, was known. Recently, 
however, evidence has been accumulating that the heavy-fermion compound UPd2Al3 (T/v = 14. 3K, T c = 1.8K, see 
Rcf. 11) is the first example where such a mechanism is responsible for superconductivity. Since this mechanism 
is pair breaking in the s-wave channel (as was already known from the impurity models mentioned above), the 
superconducting gap function must change sign as a function of k (i.e., the gap must have a node) as is the case with 
spin- fluctuation- mediated superconductivity 

An early indication that UPd2Al3 is a localised-delocalised f system came from magnetic-susceptibility measure- 
ments. Grauel et al 12 measured the dc magnetic susceptibility in UPd2Al3 and argued that the anisotropy which they 
observed arose from a tetravalent configuration of the Uranium ions, i.e. U 4+ (5/ 2 ). Since then Knight-shift 13 and 
optical 14 measurements have supported the notion that UPd2Al3 contains both localised and delocalised f electrons. 

However, the most direct evidence in favour of the dual-system hypothesis comes from inelastic neutron scattering 
(INS) and tunnelling experiments. A dispersive crystal-field excitation (magnetic exciton) was observed by Mason 
and Aeppli 15 in INS experiments within the AF phase. Later higher-resolution experiments 16-18 revealed a resonance 
feature which appears in the INS spectrum upon entering the superconducting state. This result demonstrates that 
a strong interaction exists between the localised and delocalised components of the f-electron system. Pioneering 
tunnelling experiments 19 performed on UPd2Al3-A10;E-Pb tunnelling junctions allowed experimenters to view, for the 
first time, the tunnelling density of states (DOS) of a heavy-fermion superconductor. Strong-coupling features appear 
in the DOS around lmeV close to the superconducting gap energy reinforcing the view that the exchanged bosons are 
the magnetic excitons. Taken together, INS and tunnelling experiments lead Sato and co-workers 20 to the conclusion 
that superconductivity arises in this material from an effective interaction between itinerant electrons mediated by 
magnetic excitons. Using a model two-component Hamiltonian they were able to explain qualitative features of the 
INS scattering spectrum and the superconducting tunnelling spectrum. 

Subsequently, the origin of the magnetic excitons and their global dispersion, as measured in Ref. 15, was investigated 
in more detail 21 . Using an xy-type interaction, where only the a± components of the conduction-electron spin couples 
to the magnetic excitons, the effective non-retarded pair potential was derived. The gap equations were solved within 
a weak-coupling approach demonstrating that the highest T c is obtained by an odd-parity state. In this model, 
however, the structure of the pairing amplitudes in spin space is complicated, making it unsuitable for going beyond 
the non-retarded approximation. 

In this paper, therefore, we investigate an alternative model with a simplified interaction between localised and 
itinerant 5f electrons. This interaction is of the Ising type, i.e., only the a z component of the conduction-electron spin 
can scatter magnetic excitons. In this case, the gap equations naturally divide into those for equal- and opposite-spin 
pairing, in contrast to the usual 'singlet' and 'triplet' classification that arises in spin-rotation-symmetric models. 
In the present work we treat this simplified dual model for UPd2Al 3 in a more sophisticated strong-coupling ap- 
proach using a mapping to an electron-boson Hamiltonian and solving the Eliashbcrg equations for the frequency- 
and momentum-dependent self-energy and gap functions. The level scheme which we use here, and which gives rise to 
the Ising-type interaction, is quite realistic: recent band-structure calculations based on this level scheme account well 
for experimental de Haas- van Alphen frequencies in UPd 2 Al 3 7 . We demonstrate that our model can yield a super- 
conducting transition temperature and a mass renormalisation which are consistent with experiment for reasonable 
values of the coupling constant. 

II. MODEL 

Band structure calculations based on the supposition that two of the three Uranium 5f-electrons are localised 
reproduce the observed de Haas-van Alphen frequencies in UPd2Al 3 very well 7 . These calculations suggest a level 
scheme for the localised U 5f states which we adopt here. According to the jj-coupling scheme the U 4+ (5/ 2 ) ions have 
total angular momentum J = 4. The twofold degeneracy of the ionic ground state is lifted by a crystalline electric 
field. We consider only the excitation between the non-degenerate ground state ^4) = | (| J z = 3) — \J Z = —3)) and 
the first excited state |r 3 > = \ (| J z = 3) + | J z = -3)). The CEF energy splitting is of the order A ~ 6 meV, as 
obtained from fitting the magnetic-exciton dispersion obtained from INS 15 measurements to theoretical results 21 . 

Although such a model can give rise to induced magnetism via mixing of the states |r 3 ) and |r 4 ), we take the 



view here that the underlying antiferromagnetic (AF) order is not an important consideration in the description of 
heavy- fermion superconductivity in UPd 2 Al 3 . AF order mainly leads to a folding down of the conduction bands into 
the AF Brillouin zone. Also, the magnetic-exciton dispersion in the ordered phase is not appreciably different from 
that in the paramagnetic phase, due to the fact that Tjv < A 21 . Furthermore, in mean-field theory, the AF order 
parameter will only slightly modify the superconducting T c due to a reconstruction of the conduction-electron states 
close to the AF Bragg planes. 

Accordingly we consider a three-dimensional lattice of localised 5P CEF states and itinerant f-electrons. In the 
subspace {^3), ^4)} we may write the CEF part of the Hamiltonian as 

Hqef — A S lz , (1) 

i 

where S denotes a pseudospin (S=|). In this representation we interpret the CEF ground state as having S z = -| 
and energy —A/2 and the excited CEF state as having S z = | and energy A/2. In the pseudospin representation, the 
only non-zero component of the physical total angular momentum J is the J z -component: 

Jz = l( 1 =lS x , (2) 
where 7 = 6. Then the full two-component Hamiltonian may be written as 

H = E ^ c L/k, + aJ2s iz -jJ2s i x si +s -iJ2 °i*s ix . (3) 

k.<T i iS i 

The third term in H is a nearest-neighbour superexchange interaction between localised 5f states. The last term is 
the exchange interaction between localised 5f CEF states and the delocalised 5f conduction electrons; the conduction- 
electron spin operator is 

CTi = V4( r *)°"a/#/3( r i)- ( 4 ) 

After carrying out a Holstein-Primakoff transformation (valid at low temperatures T <C A) and an additional Bogoli- 
ubov transformation involving the resulting Holstein-Primakoff bosons (see the Appendix) , H takes the form 

^E^i^ + E^K^ + VS)-/ / dv^ a {v)al^)^) (5) 

k,<7 q J 

where 



<H r ) = -j= E \ X i ( a i + a -q) 



A 



A q = -. (6) 



cj q is the dispersion of those bosons which have creation and destruction operators a q and a q respectively (see the 
Appendix). 

We define the normal electron and exciton Green's functions as follows: 

G a0 (r-T',T-r') = -(T^r^rV)) (7) 
D{v-v',t-t') = -<7Wr,TMrV)>. (8) 

Note that <^(r) ~ S x (r), as can be seen immediately by comparing the terms ~ I in Eqs. (3) and (5), and so D is 
essentially a pseudospin susceptibility. Furthermore, as in the theory of phonon-mediated superconductivity, <}f)(r) is 
real, 4>(r) = <^(r), and commutes with itself, [<f>(r), <j>(r')] = 0, conditions which allow us to use Wick's Theorem. The 
Eliashberg equations for the conduction-electron self-energy E(p,iw„) and gap function $(p,iw n ) which follow from 
the Hamiltonian Eq. (5) are 

S(p,iw„) = —^2K(p-p',iuj n -i(j' n )G(p',iuj' n ) (9) 
$(p,iw„) = P ^ ^X( P -p',^„-z^)|G(p',z^)| 2 $(p',0. (10) 
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N is the total number of lattice sites. The kernel K is given by 



We have decided not to renormalise D° by the interaction / between magnetic excitons and electrons. This interaction 
is, however, already included in D° to the extent that we model the exciton dispersion uj q by the true experimental 
dispersion. A particularly noteworthy feature of this formulation is that the effective interaction is dominated by 
its static part and is strongly peaked in q-space at the antiferromagnetic wavevector Q = (0,0,^). This strong 
dependence on wavevector contrasts with the practically wavevector-independent interaction in the usual phonon- 
mediated superconductivity. The difference arises because, in the phonon problem, the quantity corresponding to A q 
has an additional factor, w q , leading to a wavevector-independent static phonon propagator. A strong interaction 
between the collective modes of the localised moments and the heavy conduction electrons at Q = (0, 0, ^) has actually 
been observed 17 ' 18 in UPd 2 Al 3 . The electron Green's function is related to the electron self-energy via the Dyson 
equation 

G _1 (p, iu n ) = iuj n - (e p - fj,) - S(p,iw„). (12) 

The prefactor p in Eq. (10) is the expectation value of the Ising spin-spin interaction a z a z in the spin part of the 
pair wave function |x) = \S, S z ) (S here should not be confused with the pseudospin introduced earlier): 

P=(X\& Z Z Z \X)- (13) 

In the opposite-spin pair (OSP) states, 

\x) = i¥ U) - ]m (14) 

p reduces to 

p = (n\a z a z \u) = -i- (i5) 

On the other hand, for each of the equal-spin pair (ESP) states, \x) = |Tt) an d |||), p is +1. Note that spin rotational 
symmetry is broken in a maximal (Ising-type) way in Eq. (10). Consequently the usual classification of pairing states 
into an (even-parity) singlet with S = and an (odd-parity) triplet with S = 1 is no longer valid. In our model, the 
|1, 0) pair state no longer has the same energy as that of the |1, ±1) pair states. The result is a different classification 
into an ESP doublet and two OSP singlet states belonging to irreducible representations of the T>eh group (Table I). 
Finally, the conduction-band filling, defined as the ratio of the number of electrons to the number of lattice sites, is 

k 

where n k is the quasiparticle occupation factor: 

n k = n kT +n k ; = 2T^G(k,iw„). (17) 

Experimentally, the magnetic exciton has its strongest dispersion w q along the c-direction. Neglecting the weaker 
dispersion in the a-b plane allows us to reduce the three-dimensional problem to one dimension: 



E(p„iw n ) = K(p z -p' z ,i(j n -iu' n ) f ^E±-G(p' ± ,p' z ,iuj' n ) 

$(p z ,iui n ) = p^- K(p z -p' z ,iuj n -i(j' n )<f>(p' z ,iuj' n ) J ^-L\G(p' ± , P ' z ,iuj' n )\ 



(18) 
(19) 



" - ^rYA I ^n k , (20) 
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provided we can carry out the integrals along the perpendicular direction analytically. Here, N z is the number of 
lattice sites in the z-direction and A is the area of the (hexagonal) unit cell in the plane. We model the exciton 
dispersion by 



uj(q z ) = u> ex [l + (3 cos(q z )}; < /? ~ 1, 



(21) 



This form of the dispersion along the c-direction describes qualitatively the observed excitation branch 15 . Here 
2Pojex — 8 meV is the overall dispersion width and (l-(3)co ex — 1 meV is the exciton gap at the AF wave vector 
Q=(0,0,-). The quantity uj ex is a characteristic exciton energy. We choose (3 — 0.8 and ut ex ~ 5 meV ~ 60 K. The 
detailed RPA theory of the magnetic-exciton dispersion and a fit to the INS data of Rcf. 15 was given in Ref. 21. 

That sheet of the Fermi surface of UPd2Al3 thought to be most important in bringing about heavy-fermion 
behaviour 21 has the shape of a corrugated cylinder with its axis aligned along the c-direction of the hexagonal 
lattice. We therefore choose to model the electron dispersion as the sum of a strongly dispersive part in the plane ei^ 
and a weakly dispersive part in the c-direction ek z . We approximate the hexagonal unit cell in the plane (k Xl k v ) by 
a circle with radius k , chosen so that the hexagon and circle have the same area. (The areas must be the same in 
order that the maximum value of the band filling is two in both cases.) Then, assuming a parabolic dispersion in the 
plane = e±w 2 (0 < w = k±/ko < 1), we may carry out the integrals in Eqs. (18), (19) and (20) analytically. In 
the following we measure all energies, temperatures and frequencies in units of e±. The reduced Eliashberg equations 
are 



T,(p z ,iu n ) = — ^ K(p z -p' z ,iu n -iuj' n )G z (p' z ,iu' n ) 
\(T)$(p z ,iw n ) = p— ^ K(p z -p' z ,iw n -iu'JM z (p' z ,iw'J<f>(p' z ,iuj' n ) 



n 



(*' - l) 2 + (z"f 



{z'f + (z") 2 



where 



K(q z ,iv n ) 

G'^p' z ,iu> n ) 
M z (p' zl iuj n ) 



(WqJ(J ex ) + (is n /uj ex ) 2 



-I in 

2 



{z' - l) 2 + {z"f 



tan 



tan 



{z'f + {z"Y 
1-z' 



tan 



— tan 



and 



A(T) = 1 for T = T c . 



(22) 
(23) 
(24) 

(25) 
(26) 
(27) 
(28) 

(29) 



We have lumped together numerical prcfactors and the original coupling constant I into a new coupling constant g 
having dimensions of energy: 



9 = 



I 2 A ( 1 k 2 



1 



2c 2tt / u> 2 ' 



(30) 



where c is the lattice constant in the z-direction. The value of g is not calculated; rather it is considered a model 
parameter motivated by experiment. The complex number z has a real part 



and an imaginary part 



z' = - (e p ' z - fi + E'(p' z ,iw n )) 
z" =w„- Z"(p' z ,iu) n ). 



(31) 



(32) 
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TABLE I: Spin and orbital structure of the possible gap functions which are solutions of the Eliashberg equations for the dual 
model of UPd2Al3 used in this paper. 



p 


Ix) = \s, s z ) 


D 6h repres. 


spin pairing 




-1 


|o,o> = ^(ITD-UT)) 


rf (Ai 9 ) 


OSP 


cos(cp z ) 


-1 


|i,o> = ^(|U> + Ut» 


r- (Ai u ) 


OSP 


sin(cp z ) 


+1 


|i,±i> = |TT),IU> 


IT (Ai„) 


ESP 


sin(2cp z ) 



We choose a simple tight-binding form for 

e K = f cos(^); a < 1 (33) 

where a determines the degree of corrugation of the FS cylinder along c. 

The momentum convolutions in Eqs. (22) and (23) were evaluated with the aid of a fast-Fourier-transform algorithm 
on a 32 x 32 lattice. The corresponding frequency sums were carried out using the renormalisation group technique of 
Pao and Bickers 22 which allows a considerable reduction of the computational effort. Between 240 and 480 Matsubara 
frequencies were kept at each stage of the renormalisation group procedure. The renormalisation procedure was started 
at a temperature To = O.Ole^ and the frequency sum cut-off used was f2 c w 15ej_. The renormalisation procedure 
restricts us to discrete temperatures so that the point at which the condition in Eq. (29) is met must be determined 
by interpolation. The discrete temperatures were sufficiently close that a linear interpolation was adequate. 



III. RESULTS AND DISCUSSION 



Scattering by isolated impurities with CEF excitations is usually pair breaking in the s-wave channel 8 because the 
dipolar exchange interactions which break the singlet state are stronger than spin-conserving quadrupolar interactions 
between conduction electrons and CEF states. The s-wave state (where &(p z ) is nodeless) is also not favourable in 
our case of a periodic lattice of partly localised 5f-electrons which has dispersive CEF excitations. One can most 
clearly see this in the following way. In the singlet channel p = — 1. Because the kernel K(q,iv n ) is strongly peaked 
at q z = £ and v n = 0, the gap equation may then be crudely approximated in the following way: 

$(p Z) inT) ~ - C(p z , ^) <5>{p z - ^,*7iT), (34) 

where C(p z , ^ ) is a positive number. Therefore a finite gap &(p z ,inT) must change sign on translation through 
Q z = —, i.e., the gap must possess the symmetry $(p z — -) = — &(p z ), which excludes an s-wave gap function. 

However a non-s-wave superconducting state with nodes may take advantage of the exciton dispersion and become 
stable. From numerical solutions of the linearised gap equations for our model we find that the instability with the 
highest T c is accidentally doubly degenerate. The corresponding gap functions transform as the even- and odd-parity 
OSP states, cos(cp z ) and sm(cp z ). The former is the usual singlet state; the latter is the S z =0 part of the triplet 
which would appear in a theory with full spin rotational symmetry. We also find an instability in the ESP channel 
but at a much lower temperature. The corresponding pairing symmetry is &(p z ) ~ sin(2cp z ). We summarise our 
results in Table I. 

The symmetries of the orbital pair wavefunctions can be understood in the following way. By comparing our self- 
energy equations with those obtained from a four-fermion interaction we see that the effective interaction between 
quasiparticles in our theory is 

(7274|^|7i73)=/ 2 ^°< 72( t; 74 , (35) 

or equivalcntly 

V = I 2 D°a z a z . (36) 

Hence the pairing interaction in some \x) channel is 

V x (q,iis n ) = (x\a z ^\x)I 2 D (q,iu n ). (37) 

We note that our equations go beyond previous calculations 6,21 in that our formulation includes the full momentum 
and frequency dependence of the effective interaction V x (q, iis n ). 
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In the OSP channel the static interaction is 



V(q z ) = -I 2 D°(q z ) = 



PA/2 



2 ' 



(38) 



which is strongly peaked at — . The corresponding interaction in real space is therefore attractive when the quasiparti- 
cles are seperated by a lattice spacing in the c-direction. Now cos(cp z ) and sin(cpz) are peaked with equal amplitude 
in real space at z = c. Hence they are equally well suited to take advantage of the attractive part of the interaction 
and have equal superconducting transition temperatures. The ESP interaction is the negative of the OSP interaction 
and so it first becomes attractive at z = 2c in real space. The only odd-parity wavefunction which is peaked in real 
space at this position is sm(2cp z ). 

Ultimately, the degeneracy of the OSP states is a peculiarity of the particular CEF level scheme we adopt for the 
localised 5f states and of the approximate form of their exciton dispersion which we assume (viz. strongly wavevector 
dependent only in the z-direction). Both OSP states are, in fact, compatible with present experimental evidence. 
OSP states, regardless of the value of total spin angular momentum S, give rise to a reduction in the paramagnetic 
susceptibility on entering the superconducting state. This reduction comes about because oppositely paired spins 
have no magnetisation. Furthermore, such a reduction in the Pauli susceptibility leads to an upper critical field H C 2 
which, in some cases, can be smaller than that due to the Meissner effect. This effect is called Pauli limiting 23 . 
A pronounced Knight-shift reduction at Pd sites below T c (Refs. 13,24) and H C 2 Pauli-limiting 25 have both been 
observed in UPd 2 Al 3 . Note also that both OSP gap functions have node lines perpendicular to the c-direction. In the 
even-parity case they are located at k z = ±^Q Z , i.e., at the Bragg planes and zone boundaries of the AF Brillouin 
zone, whereas in the odd-parity case they are located at k z —0. Experimentally, the existence of node lines was inferred 
from NMR relaxation-rate 26 ' 27 , specific-heat 28 and thermal-conductivity measurements 29 . These experiments did not, 
however, locate the node position along the c-direction and so the correct gap symmetry in UPd2Al 3 presently remains 
an open problem. To distinguish between the possible nodal gap functions, it is very important to perform field-angle 
dependent measurements of the specific heat or of the thermal conductivity at low temperature. As proposed in 
Ref. 30, such measurements may be able to locate the node line in q-space. 

The results of our numerical calculations of the superconducting transition temperature in the ESP and OSP 
channels are shown in Fig. 1. The superconducting instability occurs first in the OSP channel. We have carried out 
the calculations with u ex = 0.01ej_. This exciton energy corresponds to about 60K for an electronic bandwidth of 0.5 
cV. One therefore obtains a transition temperature in the OSP channel close to the experimental value 11 (T c = 1.8K) 
with a dimensionless coupling constant g/e± of about 0.5. At this point we remind the reader that this temperature 
is in the range of validity of our theory T c -C u) ex ~ A. 

That contribution to the mass renormalisation m* /mj (mi is the band mass) arising from the momentum depen- 
dence of the real part of the self energy is small and so m* /nib simply reduces to the Eliashberg renormalisation factor, 
Z(p,iu> n ), which is practically momentum-independent. The weak dependence of Z(p,iu> n ) on p follows from the 
combined effect of a strong exciton dispersion along the c-axis and a weak electron dispersion in the same direction. 
Our results for the mass renormalisation are shown in Fig. 2. Note that the mass renormalisation is approximately 
linear in g ~ I 2 (see Ref. 31). Ah initio calculations of the band masses nib/mo (mo is the free electron mass) for the 
7 ring and [3 ring of the cylindrical Fermi surface have been carried out by confining two of the three 5f electrons to 
the Uranium ions 7 . These values, in conjunction with a mass enhancement of m'/mt of about 10, yielded values of 
m* /mo which were in good agreement with experiment 32 . We find that such a large value of the mass renormalisation 
corresponds to a value of the coupling constant (g/e±_ ~ 2) which is roughly a factor of 4 larger than that required to 
reproduce the transition temperature (g/e± sa 0.5). One should note that the AF long-range order may act to reduce 
T c , an effect which has not been considered here. Therefore the proper g may indeed be larger than 0.5. 

We have also carried out our calculations in the case that the Fermi surface is purely cylindrical (a = 0). The 
results are practically unchanged. This case, however, has the appealing feature that the bracketed factor in the 
definition of the coupling parameter (see Eq. 30) is related in the following simple way to the (constant) density of 
states per spin, N(u), 



g = 




(39) 



The dimensionless coupling constant can now be written 



1 2 
9/e± = g ( A i A a) A 3, 



(40) 



co = O.Ole,; a = -0.1; (3 = 0.8 
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FIG. 1: The dependence of the superconducting transition temperature T c on the coupling constant g [Eq. (30)]. The band 
filling n is 0.6. 
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FIG. 2: The dependence of the mass renormalisation m* /m& on the coupling constant g [Eq. (30)]. The band filling n is 0.6 
as in Fig. 1. m* min and mj na;r are the extremal effective masses over all values of q z . vs\* min occurs at q z — 0; vs\* max occurs 
at q z = ^. m* is practically constant over the Fermi surface for the range of values of g which we consider in this work 
(0.5 < g/e ± < 4). 
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where the dimensionless constants Aj are defined by 

Ai = IN{ijl) (41) 
A 2 = A/u ex (42) 

A3 = Amry ™ 

Of these three parameters, the first, Ai, is not presently known for the Uranium 5f systems. The second, A 2 , is about 
unity. The third, A3, is approximately 167, if we make the estimate N(n) ss leV -1 per cell per spin. The value of the 
coupling parameter we require to get the experimentally observed transition temperature then leads to the estimate 
Ai s=a 0.08. A value of Ai which is a factor of 2 larger would yield a mass renormalisation in good agreement with 
experiment, but then T c would be larger than the experimental value (1.8K) by roughly a factor of 1.6. As already 
mentioned, this discrepancy may be alleviated by the effects of AF order. The values for Ai obtained here for a dual 
5f compound are of the same order of magnitude as for 4f systems 9,31 . 



IV. CONCLUSIONS AND OUTLOOK 



Motivated by experiments on UPd2Al3, we have examined, for the first time, a magnetic-exciton model of supercon- 
ductivity within a strong-coupling retarded framework. The model naturally explains the strong interaction between 
the collective modes of the localised moments and the heavy conduction electrons which is observed at Q = (0, 0, -) 

in UPd 2 Al 3 . Solutions of the Eliashberg equations show that the model favours a superconducting instability for the 
even- or odd-parity OSP states. Each state has line nodes perpendicular to the c-direction and total z-component of 
spin S z — 0, characteristics which are compatible with measurements of UPd2Al3. We find that a superconducting 
transition temperature and mass renormalisation, each in good agreement with their experimental values in UPd2Al3, 
can be obtained using reasonable values of the parameters in the theory. Taken together, these results strengthen the 
argument in favour of a magnetic-exciton system in UPd2Al 3 . 

The evidence for dual f-systems in other Uranium compounds does not weigh as heavily as it does in the case of 
UPd2Al3. In UPt3 the number of itinerant f-clcctrons remains controversial. However band-structure calculations 
based on the assumption that two of the three 5f-electrons are localised reproduce the observed de Haas- van Alphen 
frequencies as well as the anisotropic heavy electron masses in this compound very well 7 . 

Recently unconventional superconductivity with split-T c superconductivity has been discovered in the Pr-based 
skutterudite PrOs4Sbi 2 33,34 . Thermodynamic and transport measurements suggest that the observed heavy-fermion 
state arises from the interaction of electric quadrupole moments of the CEF-split Af 2 states of Pr 3+ with the conduc- 
tion electrons 33 . The ground state is probably a singlet and the first excited state is probably a triplet, approximately 
0.5 meV higher, with large off-diagonal quadrupolar matrix elements connecting the two states. The largeness of these 
matrix elements, together with the very small singlet-triplet splitting, strongly suggests an effective mass renormalisa- 
tion via virtual CEF excitations. Finally, experiments reveal that the superconducting order parameter is anisotropic, 
making this compound another candidate for CEF-exciton-mediated superconductivity, but this time of quadrupolar 
nature. 

UPd 2 Al 3 provides a strong motivation for the study of dual f-electron systems such as the one discussed in this work. 
We have considered a dual model in which all parameters are fixed except the coupling constant. It will be interesting 
to investigate the dependence of the transition temperature on conduction-band filling and on the form of the exciton 
dispersion, and to look for a mechanism that can discriminate between the even- and odd-parity superconducting 
OSP states. 
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APPENDIX 

In this appendix we sketch the derivation of Eq. (5). Our starting point is the set of Holstein-Primakoff transfor- 
mations (valid for T<A), 



S x = ~+aU (A.l) 



2 

S + -Z a a t] => S x = J(a + at) (A.2) 

Inserting these expressions into the second and third sums in Eq. (3) and Fourier transforming, we arrive at the 
following form, 



^l(q)«q a q _ ^-y^(«q a -q + a q a -q) ~ 4 U *<*Zfi C k+<i0 K + fl -q)> ( A ' 3 ) 

q L J kq 



where 

Wi(q) = A-W 2 (q). (A.4) 

The first sum can be diagonalised 35 , while at the same time introducing only a factor A q in the second sum. To do 
this we use the Bogoliubov transformation 



(A.5) 



where u k = cosh(#t) and v k — sinh(#k)- If #k is chosen such that 

-W 2 (k)/Wi(k) = tanh(20 k ), (A.6) 
then (ignoring constant additive terms) the first sum in Eq. (A. 3) takes the form, 

^> q (a q a q + l/2) (A.7) 



where the dispersion w q satisfies 
Finally, it is easy to show that 
where 



"1 = Wi(q)-W|(q). (A.8) 

fl q + a -q = \( a q + a -q), (A. 9) 



A q = cosh((9 q ) - sinh(0 q ). (A.10) 
Using Eqs. (A.4) and (A.6), together with the identity 



A ' 

we can then show that 



A-2W 2 (q) = ^, (A.ll) 



A q = -. (A.12) 



OJ, 



q 



The interaction term in Eq. (5) follows immediately. 
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